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In Ref. [1] the compression-relaxation mechanism 
(CRM) for a graphene sheet (GS) was reported to be irre- 
versible and resulted in static ripples on GS, such that for 
T <Tc the free-energy of the rippled GS is smaller than 
that of roughened GS. We will point out several techni- 
cal difficulties, such as the use of the relation AA = (W) 
for the free energy calculations and the definition of the 
rough state, with their simulations [2J. We show that 
(at T <Tc) their introduced rough state suffer boundary 
stress, thus is a rippled state and their obtained inequal- 
ity A r ippi e <A roug h is no longer valid. Therefore the intro- 
duced mechanism is reversible. Furthermore, from theo- 
retical point of view for an infinitely slow rate of CRM Q , 
the relation AA = (W) (in common non-equilibrium sim- 
ulations) is allowed Q. In this case authors of Ref. [f[ 
must report the used infinitely small rate and justify how 
it is in practice valid. We show AA j£ (W) for the system 
that was used in Ref. 

To define the rippled state it was stated in Ref. [l|: 
'we might state that any primary stress on GS, for ex- 
ample, in its preparation in the experiments, can con- 
struct ripples that will survive during the experimental 
measurements' . According to this statement, we assume 
that each GS (nanoribbons) with boundary stress is in 
the rippled state (regardless of the shape of the ripples 
and the other corresponding properties). 

We simulated the same sample that was used in 
Ref. [l|. A GS with 80x40 atoms by using the same 
method at T=55 K, the strain rate for CRM is 0.0117/ns. 
The time step is 0.5 fs and after 25 ps (for thermal equi- 
librium, i.e. state '/'), we compressed the system for 325 
ps and then relaxed it by moving back to its initial posi- 
tion (state 'F') for other 325 ps. Obviously this is done 
by applying forces on the two rows of atoms at the (lon- 
gitudinal) ends. Fig. [Ua) shows 8 trajectories which are 
the variation of the work done on the GS versus time (8 
simulations were done with different initial conditions) 
and the variation of ((h 2 )} (Fig.QJb)). The dashed curve 
is (W) which its variation is related to the dissipated 
work [4|. As we see (from Fig. [IJb))((/i 2 )) (after relax- 
ation) decreases with time. Therefore the amplitude of 
the final ripples decreases after CRM. Fig. [TJc) shows 
the change in the free energy using Jarzynski equal- 
ity (i.e. AA = -^- 1 ln(exp(-/3I^))), (W) and AA - 
(W)(= -0.5f3al) and clearly AA ^ (W). The difference 
AA — (W) is just the dissipated work which is associated 
with the increase of entropy [!, H| . Practically, the latter 
difference must be considered even by very small rate of 
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FIG. 1: (Color online) (a) Variation of the total works for 
8 different trajectories (solid curves) and the average work 
(dashed curve), (b) ((h 2 )} vs time (inset is zoomed region 
after relaxation), (c) free energy change, the average work 
and dissipated work and (d) boundary stress during CRM. 



CRM which is in contrast to the assumption AA = (W) 
in Ref. [I]. 

The size of the system at state T is assumed to be 
equal to a flat GS (it is due to performing NVT ensem- 
ble in Ref. [I]) for a sample size with 80 x 40 atoms 
('rough state' in Ref. [I|). Fig. [IJd) shows the varia- 
tion of (absolute value of) boundary force per width (i.e. 
stress) versus CRM time (the averaged value over 8 sim- 
ulations). As we see from the top panel in Fig. [IJd) at 
state ('rough state' in Ref. [I]) the system is not free 
from stress and it feels pre-stress around 4 GPa which is 
not negligible @ (the reason is that the GS does not have 
the optimum size @ at '/'). Here we emphasize that the 
behavior of boundary stress around points (and 'F') is 
always similar to that one we presented here. After some 
compression steps the free energy curve shows the first 
minimum, i.e. state 'O'. At this state the system is al- 
most free of boundary stress (Fig. [IJd)). Moving back to 
the initial position by stretching the system, after passing 
another minimum [7| we obtain the state l F', which is not 
free from boundary stress too ('ripple state' in Ref. [I|). 
Therefore both the initial and final state in CRM are 
rippled states @. Moreover, the state 'O' should be con- 
sidered as the true 'rough state' in Ref. [I], where there 
is no boundary stress and the relation Ao — Ap < 
holds @. The latter inequality implies that the state 
'O' is much stable than the rippled states 'F' and 
In Ref. [I] the inequality Ap — Ai is referred to as the 
difference between rippled (with stress) and the rough 
state (without stress) i.e. difference between two dashed 
horizontal lines in Fig. [IJc), however we see that this is 
only the difference between two rippled states, i.e. 
and l F' . The final remark is that in Ref. [I], there is no 
physical reason given for the irreversibility except that 
Ap — Ai < 0, which we have shown this is not valid. 
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Note that the existence of the boundary stresses at the 
states '/' and 'F' are independent of other parameters 
(e.g. rate for CRM, amount of compression, size of the 
system and temperature Q). 

In summary, the true 'rough state' has lower free en- 



ergy with respect to the rippled state. Using the Jarzyn- 
ski equality is necessary when calculating the free energy 
change in the CRM of GS even for very slow evolution. 
Both initial and final states in the simulations of Ref. [1] 
are rippled states thus the CRM is reversible. 
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